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The local atomic structure of an amorphous NiZr2 alloy was investigated using the anomalous 
wide-angle x-ray scattering (AWAXS), differential anomalous scattering (DAS) and reverse Monte 
Carlo (RMC) simulations techniques. The AWAXS measurements were performed at eight different 
incident photon energies, including some close to the Ni and Zr K edges. From the measurements 
eight total structure factor S{K, E) were derived. Using the AWAXS data four differential structure 
factors DSFi(K, Em, E„) were derived, two about the Ni and Zr edges. The partial structure 
factors 5Ni-Ni(-?s'), 5Ni-Zr(-ft') and Szr-Zr{K) were estimated by using two different methods. First, 
the S{K,E) and DSFi{K, Em, E„) factors were combined and used in a matrix inversion process. 
Second, three S{K, E) factors were used as input data in the RMC technique. The coordination 
numbers and interatomic distances for the first neighbors extracted from the partial structure factors 
obtained by these two methods show a good agreement. By using the three-dimensional structure 
derived from the RMC simulations, the bond-angle distributions were calculated and they suggest 
the presence of distorted triangular-faced polyhedral units in the amorphous NiZr2 structure. We 
have used the Warren chemical short-range order parameter to evaluate the chemical short-range 
order for the amorphous NiZr2 alloy and for the NiZr2 compound. The calculated values show 
that the chemical short-range order found in these two materials is similar to that found in a solid 
solution. 

PACS numbers: 61.10.Eq, 61.43. Bn, 05.10.Ln, 81.15.Cd 



I. INTRODUCTION 

Amorphous materials have a great potential for appli- 
cation in technological devices, but their uses are limited 
due to several factors. One of them is the difficulty to 
obtain information about their atomic structures. Since 
most chemical and physical properties of these materi- 
als are defined by the first neighbors, knowledge of their 
structures becomes necessary. Due to the lack of transla- 
tional symmetry, the determination of their atomic struc- 
tures involves the combination of different diffraction and 
spectroscopic techniques as well as simulations and mod- 
eling. 

The structure of an amorphous binary alloy is de- 
scribed by three pair correlation functions Gij{r), which 
are the Fourier transformation of the three partial struc- 
ture factors Sij{K). The total structure factor S{K,E), 
which can be derived from scattering measuiements, is 
a weighted sum of these three Sij{K) factorad. Thus, to 
determine the three partial Sij {K) at least three indepen- 
dent S{K, E) factors are needed. Usually, the following 
methods have been used to obtain these S{K, E) factors: 

1. Isomorphous substitution: in this method one com- 
ponent of the amorphous alloy is substituted by a 
(expected) chemically similar one. For example, in 
amorphous Ni-Zr alloys the Zr atoms are replaced 



by Hf atoms or the Ni atoms are partially substi- 
tuted by Co or Fe atom^. This method is not very 
accurate since one cannot be sure of the chemical 
similarity of the two elements involved. 

2. Isotope substitution: in this method one compo- 
nent of the amorphous alloy is substituted by an 
isotope or by an isotopic mixture and the scat- 
tering measurements are carried out using neutron 
sources. For example, in amorphous Ni-Zt, alloys 
the Ni atoms are replaced by a Ni isotopeQ. The 
advantages of this method are limited by the avail- 
ability of adequate isotopes. 

With the development of the synchrotron radiation 
sources, anomalous wide angle x-ray scattering (AWAXS) 
and differential anomalous scattering (DAS) techniques 
became available for structural study of multicomponent 
disordered materials. AWAXS utilizes an incident radi- 
ation that is tuned close to an atomic absorption edge 
so that it interacts resonantly with the electrons of that 
particular atom. The atomic scattering factor of each 
chemical component can therefore be varied individually 
and the chemical environment about each component in 
the material can be investigated. Thus, in the case of an 
amorphous binary alloy, by using only one sample three 
independent S{K, E) factors can be obtained. However, 
the matrix formed by the weights of the three S{K, E) 
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factors is ill-conditioned, becoming difficult the determi- 
nation Qfrjhe three Sij{K) factors. 

FuossQi3 tried to overcome this difficulty by imple- 
menting |-|tlie DAS approach, which was proposed by 
SchevickOlD. The DAS approach consists in taking the 
difference between the scattering patterns measured at 
two incident photon energies just below the edge of a 
particular atom, so that all correlation not involving this 
atom subtract out since only the atomic scattering factor 
of this atom changes appreciably. Later ,JI)e Lima et al.u, 
following a suggestion made by MunroLl, combined the 
differential scattering factors DSFi{K, Em, En) and the 
S{K, E) factors. They observed that this combination re- 
duces the conditioning number of the matrix formed by 
the weights of these factors, allowing more stable values 
of Sij {K) to be obtained. 

•esse Monte Carlo (RMC) simulation techni- 



and 



queEjUOEfl has been successfully used for structural 
modeling of amorphous structures. S{K,E) factors 
or their Fourier transformations, the total reduced ra- 
dial distribution functions G{r), can be used as inj 
data. AjjpJications of this. <iPcl | i^n| ^ q^C || t(j 
crystalline^ and amorphoustatZll 
materials are described in the literature. 

In this paper, we report the determination of the 
three Sij{K) factors obtained for amorphous NiZr2 al- 
loy by two independent ways: by making a combination 
of AWAXS and DAS techniques and also by combining 
AWAXS and RMC simulations techniques. 



{f{K,E)) = Y.cJ^{K,E), 

i 

{f{K,E)f - \^cJ,{K,E) 



Here, Ci is the concentration and fi(K,E) ~ fo{K) + 
f {E) + if" (E) is the atomic scattering factor of the com- 
ponent of type i and f'{E) and f"{E) are the anomalous 
dispersion terms. The total reduced distribution function 
G{r) is related to the S{K, E) factor through the Fourier 
transformation 
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G{r) = - J K [S{K, E) - l] sin(i^r) dK , 



while the partial reduced distribution functions Gij{r) 
are related to the Sij (K) by means of the Fourier trans- 
formation 
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G^j ir) = - K [S^j (K) - 1] siniKr) dK . 
^ Jo 



II. THEORETICAL BACKGROUND 

A. Total and partial structure factors 

According to Faber and ZimanEI, the S (K, E) factor is 
obtained from the scattered intensity per atom la {K, E) 
as follows: 



From Gij (r) the partial radial distribution functions 
RDFij(r) can be calculated by 



RDFy(r) = Airpocjr^ + rGy(r) . 



SiK,E) = 



Ia{K,E)^[{P{K,E))-{f{K,E)f] 
{f{K,E))^ 

n 71 
i=l 3 = 1 



(1) 



where K is the transferred momentum, E is the incident 
photon energy, Sij{K) are the partial structure factors 
and Wij {K, E) are given by 



{f{K,E)Y 



Interatomic distances are obtained from the maxima of 
the Gij (r) functions and coordination numbers are calcu- 
lated by integrating the peaks of the RDF^ (r) functions. 



B. Differential structure factors and differential 
radial distribution functions 

By using the Faber and Ziman formalismB, the 
DSFi{K, Em, En) factor around the component i is ob- 
tained as follows 
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DSF,{K,E^,En) = 



[la{K, Em) - IgjK, En)] - [LSjK, E„,) - LSjK, E^)] 



{f{K,E„,))^-{f{K,E,,))^ 
= Y,U^JiK,E,r,,E,,)S,,{K), 



(3) 



where 



is the Laue scattering term and 



LS{K,E) ^ {f{K,E)) - {f{K,E)f 



U^j{K,Era,En) 



c,cj [MK, E„,)fAK, E,n) - MK, En)fj{K, E^)] 
{f{K,Em))'^{fiK,E,,))^ 

I 



(4) 



For i ^ j the weight Uij{K, Ej^, En) must be multiphed 
by 2. The differential distribution function DDFi{r) can 



be calculated by 
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DDFi{r) ^ iirpr^ + — K[DSFi{K, £;„,£;„)- l] sin(is:r) dK . 
Jo 



r 



For a binary amorphous alloy, as compared to the to- 
tal radial distribution function RDF(r), which is a sum 
of three RDFij(r) functions, this function is the sum of 
only two RDFy (r) functions and is sensitive only to the 
environment around the atoms of type i. In this way it 
is similar to EXAFS, though there are important diffec.- 
ences due to the different sections of iC-space measuredn3. 



as seen in Eqs. || and ^ is a weight which allows 
us to vary the contribution of any S!Fi{K) factor. Here, 
Pi = 1, V£, and M — 12. By taking the derivatives of the 
J-{K) function with respect to the Sij(K) factors and 
putting them equal zero a linear system formed by three 
equations can be found. These three linear equations can 
be written as a matrix equation through 



C. The matrix inversion method 

In order to determine the three Sij{K) factors using 
the eight S{K,E) and four DSFi{K,Em,En) factors, a 
least-squares method was used. First of all, we define the 
function 



S ^ WS: 



and by inverting this matrix equation the Sij {K) can be 
found. There are several methods to invert a matrix, 
and we have adnpted the Singular Value Decomposition 
(SVD) methodET 



Mr- 2 2 

For 1 < ^ ^ 8, STi{K) indicates the S{K,Ei) factors 
and for 9 < £ ^ 12 it stands for the DSFf{K,E„,,En) 
factors, as defined by Eqs. |l| and |^. In the same way, for 
1 < ^ ^ 8, h\^{K) indicates the Wij{K,Ei) weights and 
for 9 < £ 12 it stands for the U^j{K, E„,, En) weights. 



D. The RMC method 



The basic idea and the algorithm 



.standard 
and its 
reported in the 
In the RMC 



RMC method are described elsewher 
applicaticm^ba-^dijBEejfe^ 
literature0mllZlffl@BBr 

procedure, a three-dimensional arrangement of atoms 
with the same density and chemical composition of the 
alloy is placed into a cell (usually cubic) with periodic 
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boundary conditions and the G'^'^^(r) functions corre- 
sponding to it are directly calculated through 



'iirpor'^Ar ' 

where rt^^^(r) is the number of atoms at a distance 
between r and r + Ar from the central atom, averaged 
over all atoms. By allowing the atoms to move (one at 
each time) inside the cell, the G^'^'^ (r) functions can be 
changed and, as a consequence, the Sij (K) and ^i^^c (-^^ 
factors are changed. Thus, the S^^'^{K,E) factor is 
compared to the S{K, E) factor in order to minimize the 
differences between them. The function to be minimized 
is 

i—l 

where the sum is over m experimental points and 5 is 
related to the experimental error in S{K,E). If the 
movement decreases ip'^, it is always accepted. If it in- 
creases it is accepted with a probability given by 
exp(— A'0^/2); otherwise it is rejected. As this process 
is iterated '0^ decreases until it reaches an equilibrium 
value. Thus, the atomic configuration corresponding to 
equilibrium should be consistent with the experimental 
total structure factor within the experimental error. By 
using the G^^^(r) functions and the S^^'^{K) factors 
the coordination numbers and interatomic distances can 
be calculated. In addition, the bond-angle distributions 
can also be determined. 



III. EXPERIMENTAL PROCEDURE 
A. Sample and holder 



monochromator and a Si:Li energy sensitive detector 
which is sufficient to discriminate the large Kq, resonant 
Raman/fluorescence signal when incident photon energy 
is tuned close to the Ni and Zr K edges but insufficient 
to distinguish the small elastic, and Compton sig- 
nals. The Si (220) and Ge (400) crystal monochroma- 
tors were used to collect data around the Ni and Zr 
edges, respectively. The sample was sealed, under vac- 
uum, into a cell containing a large kapton window fixed 
around the diffractometer horizontal rotation axis to re- 
duce the air scattering. In order to take into account the 
time-decrease of the beam, the incident intensity Iq was 
measured by placing a Nal photomultiplier tube between 
the monochromator and the diffractometer. Absorption 
measurements were performed on the Ni and Zr K edges 
for energy calibration. For each incident photon energy 
listed in Table |, the energy windows for Kq, , resonant Ra- 
man/fluorescence, K^, Compton, and elastic signals were 
defined by measuring the scattered intensities at two dif- 
ferent angular positions: the beginning of the windows 
was defined at the lowest scattering angle and their ends 
at the highest one. Scattering patterns were collected us- 
ing the symmetrical reflection geometry. To improve the 
quality of the data, at least four scattering measurements 
were performed at the incident photon energies shown in 
Table |. 



TABLE L // and /" values used here. 



Energy (eV) 


fk 




/zr 


/zr 


8210 


-3.895 


0.491 


-0.347 


2.167 


8330 


-7.333 


0.658 


-0.367 


2.167 


9659 


-0.964 


3.051 


-0.586 


1.629 


11867 


-0.063 


2.158 


-0.950 


1.126 


16101 


0.254 


1.287 


-1.985 


0.648 


17398 


0.269 


1.121 


-2.999 


0.563 


17987 


0.269 


1.057 


-7.371 


0.734 


19200 


0.264 


0.940 


-1.953 


3.315 



Amorphous NiZr2 thin films were prepared on a NaCl 
cooled substrate by the sputtering triode technique at 
the Argonne National Laboratory. In order to eliminate 
the substrate a metallic holder with a central rectangular 
hole was glued on the thin film, and the set was immersed 
in water to dissolve the NaCl substrate. The thickness 
and density of the thin films were 2 /xm and 7.32 g/cm^, 
respectively. 

B. Apparatus and data collection 

The AWAXS scattering experiments were performed 
some years ago on the D4 beamline at LURE (Orsay - 
France) using the DCI synchrotron radiation source. All 
measurements were performed using a two-circles diffcaG-, 
tometer in the vertical plane described elsewheretJ'Ell 
equipped with a two crystal Si (220) or Ge (400) 



C. Determination of the anomalous scattering 
factors 

In order to interpret the scattering data correctly, the 
real and the imaginary parts /' and /" of the atomic 
scattering factor were determined a<jejirately following a 
procedure described by Dreier et allS3 X-ray absorption 
coefficients were measured over a broad energy range near 
the Ni and Zr K edges on the sample and /" was calcu- 
lated using the optical theorem. The absorption measure- 
ments were performed on the XAS beamline at LURE, 
in transmission mode, and two standard ionization cham- 
bers were used as the detecting systeni,-. Outside the re- 
gion of measurement, theoretical valuesES of /" were used 
to extend the experimental data set over a larger energy 
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range and /' was calculated using the Kramers-Kronig 
relation. For the measurements away from the K edges, 
the /' ard-/" values were taken from a table compiled 
by SasakiEJ. The obtained values for the incident photon 
energies 8330 eV and 17987 eV are listed in Table | to- 
gether with Sasaki values. The atomic scattering factor 
away from the K edge, fo{K), of neutral Ni and Zr atoms 
were calculating accocding to the analytic function given 
by Cromer and Manr£j. 



D. Data analysis 

Each scattering pattern was treated separately. Due to 
the diffractometer characteristics, the polarization cor- 
rection was disregarded. The measured scattered inten- 
sities away from the Ni and Zr K edges were corrected 
for detector non-linearities and reabsorption effects be- 
ing then put on a per atom scale and the Compton scat- 
tering eliminated. For those measured close to the Ni 
and Zr edges, besides the corrections above, the con- 
tribution of the Kjs fluorescence to the elastic scatter- 
ing was eliminated. This contribution was evaluated, 
after corrections for reabsorption effects, by using the 
K^/Kq ratio, which was measured at an incident pho- 
ton energy above the edge, because they were completely 
separated. Two method were tested to put the scat- 
tered intensities on a pffl .. ati Om scale: that described 
by Krogh-Moe-Norr 



high-angle procedure 



„36,g7j|3S| g^jjj^ ^]-^g Qj^g based on the 
3j3S| Pqj. ^j^g elastic scattered in- 
tensity around the Zr edge, for which the K data range 
is larger (iCmax ~ 13 A~^), both method gave the same 
result. The K data range was reduced and the Krogh- 
Moe-Norman method was then applied. The normalized 
signal was identical to the one obtained previously. So, 
we have kept the Krogh-Moe-Norman method to put the 
scattered intensities on a per atom scale including those 
around the Ni edge, for which the K data range is smaller 
(^max ~ 8 A~^). The good quality of the total structure 
factors shown in Fig. |l| and ^ corroborates the accuracy 
of the procedure used. The Compton scattering contri- 
bution was calculated according to the analytic approxi- 
mation given by Palinkagi3. 



IV. RESULTS AND DISCUSSION 

A. Total structure factors, differential structure 
factors and differential distribution functions 

Figures |^ and || show the S{K, E) factors for the inci- 
dent photon energies described in Table |. The S{K, E) 
factors for the energies greater than 16101 eV fsee Fig. ||) 
show a pre-peak located at about K — 1.65 A^^, which 
reaches the maximum intensity close to the Zr edge. It 
is absent in the S(K,E) factors for the smallest energies 
(see Fig. |l|). Since S{K, E) is a weighted sum of Sij{K), 
we calculated the weights Wij {K, E) for the iSni-ni {K) , 



^Ni-Zr (-f'') and Szr-Zr{K) factors. Their contributions for 
the S{K,E) factor at 8330 eV are about 4%, 32% and 
64% while for the S{K, E) factor at 17987 eV they are 
about 10%, 43% and 47%, respectively. Thus, the deter- 
mination of Ni-Ni first neighbors is more difficult than 
for the Ni-Zr and Zr-Zr first neighbors. 




FIG. 1: S{K,E) factors for energies Ex = 8210 eV (A), E2 = 
8330 eV (B), E3 = 9659 eV (C) and E4, = 11867 eV (D). 




FIG. 2: S{K,E) factors for energies Es = 16101 eV (E), 
Eg = 17398 eV (F), £7 = 17987 eV (G) and Es = 19200 eV 
(H). 

We have used the scattered intensities on a per atom 
scale obtained at the incident photon energy pairs (8210, 
8330 eV), (8330, 9659 eV), (16101, 17987 eV) and (17987, 
19200 eV) to calculate four DSFi{K, E,n^n) factors 
around the Ni and Zr atoms. Figures Sand Q show them 
and their DDFi{r) functions. Figure]^ shows that the 
chemical environment around Ni and Zr atoms is very 
different. For instance, it is seen a broad shoulder on 
the right side of the first peak in the DSFj^i(K, Em, En) 
factors which is absent in the DSFzr{K, E„i, En) fac- 
tors. The DDFzr{r) functions show the first coordi- 
nation shell split into two subshells, located at about 



r = 2.64 A and 3.23 A, while the DDFm{r) functions 
display the first coordination shell located at a mean dis- 
tance r = 2.77 A. The contribution of the Ni-Ni and Ni- 
Zr pairs to them is not resolved. This can be related to 
the small K data range of the measured scattered inten- 
sities. The weight Um-m{K) to the DSF^iiK, E^n, En) 
factor is about 25%, and it is almost four times larger 
than the weight W^Ni-Ni(^) for the S{K,E) factors, al- 
lowing a more realistic determination of the number of 
Ni-Ni first neighbors. 



10 F ^ 




2 4 6 8 10 12 14 

K(k') 



FIG. 3: DSFi{K,Em,En) factors. DSFm{K, Em, E„) for 
energies {Ei,E2) (A) and iE2,E3) (B). DSFzr{K, E^, E„) 
for energies (.Es.St) (C) and (E't.^s) (D). 

By considering the atomic radius of Ni and Zr atoms, 
the subshells located at the distances r = 2.64 A and 3.23 
A in the DDFzr{r) functions were directly attributed to 
Zr-Ni and Zr-Zr correlations. This fact allowed us to de- 
convolute their first shells and, consequently, also the first 
shells of the DDFj<ii{r) functions, by assuming a Gaus- 
sian shape for them. To do this, we explicitly took into 




I I I I I I I I I I I I I I I I I I I I I 

123456789 10 11 
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FIG. 4; DDFi{r) functions. DDFm{r) for energies {Ei,E2) 
(A) and {E2, E-i) (B). DDFzr{r) for energies {E5, E7) (C) and 
{Et,Es) (D). 
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TABLE II: Structural parameters determined for the a-NiZr2 
alloy. 



Dccoiivolution of DDFi{r) functions 


Bond Type 


Ni-Ni 


Ni-Zr 


Zr-Ni 






Zr-Zr 




1.3 


8.4 


4.2 






11.6 


r (A) 


2.67 


2.67 


2.76 






3.21 


Matrix inversion 


Bond Type 


Ni-Ni 


Ni-Zr 


Zr-Ni 






Zr-Zr 


N 


2.3 


7.5 


3.8 






10.8 


r (A) 


2.64 


2.77 


2.77 






3.24 


RMC 


Bond Type 


Ni-Ni 


Ni-Zr" 


Zr-Ni' 






Zr-Zr 


N 


3.2 


6.9 2.3 


3.5 1.1 






10.1 


r (A) 


2.68 


2.73 3.62 


2.73 3.62 






3.25 






Crystalline 


NiZr2 compound 




Bond Type 


Ni-Ni 


Ni-Zr 


Zr-Ni 






Zr-Zr' 


N 


2.0 


8.0 


4.0 1.0 


2.0 4.0 4.0 


r (A) 


2.63 


2.79 


2.79 2.82 


3.17 3.30 3.47 


Amorphous NissZrgs alloyH 


Bond Type 


Ni-Ni 


Ni-Zr 


Zr-Ni 






Zr-Zr 


N 


3.3 


8.6 


4.8 






11.0 


r (A) 


2.45 


2.85 


2.85 






3.30 


Amorphous Ni36Zr64 alloy 






Bond Type 


Ni-Ni 


Ni-Zr 


Zr-Ni 






Zr-Zr 


N 


2.3 


7.9 


3.9 






9.1 


r (A) 


2.66 


2.69 


2.69 






3.15 



"There are 9.2 Ni-Zr pairs at (r) = 2.95 A. 
''There are 4.6 Zr-Ni pairs at (r) = 2.95 A. 
"^There are 11 Zr-Zr pairs at (r> = 3.29 A. 



account the K dependence by Fourier transforming into 
the K space the Gaussian distributions multiplied by the 
weights Uij{K, Em, En) and back Fourier transforming 
over the same K range used to calculate the DDFi{r). 
The structural parameters extracted from the deconvo- 
lution are listed in Table ||. 

B. Partial structure factors obtained from the 
matrix inversion and RMC methods 

1, Partial Si j{K) factors obtained from the matrix 
inversion method 

By considering the eight S {K, E) and the four 
DSFi{K, Em, En) factors, a reduction of a factor of 15 
(at K = 1.65 A-i) and 7 {at K = 7.65 A"!) in the 
conditioning number of the matrix W was reached when 
compared to the use of the eight S {K, E) factors only. 
To obtain the three Sij{K) factors, the S{K,E) and 
DSFi{K, Em, En) factors were restricted to i^max = 7.65 
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A~-^. The iSNi-Ni(^), Sm-Zi{K) and Szi-ZriK) factors 
(thin Hnes) obtained are shown in Fig. ^. From this 
figure, the 5Ni-Zr(-fi') and Szr-Zr{K) factors seem to be 
of good quahty, specially between K = 1.5 and 
K = 3.5 A~^, where the conditioning number presents 
the largest changes. On the other hand, the 5Ni-Ni(-^) 
factor displays a minimum apparently without physical 
meaning in this region. The S^i.zr{K) factor shows a 
minimum at K ^ 2.2 A-i, which was already observed 
in the DSF^i{K, E^^, E^) factors, while the Szr-7,r{K) 
factor displays a shoulder at this K value. Although the 
'^Ni-Ni [K) factor is very noisy, it is possible to see a weak 
pre-peak at about if = 1.5 A~^ which was already ob- 
served in the S{K, E) factors around the Zr edge and also 
some weak halos beyond K — 3.5 A~^. 



14 . 




FIG. 5: <SNi-Ni(if), St^i-zr{K) and Szr-Zr(K) factors obtained 
from the matrix inversion (thin lines) and RMC simulations 
(thick lines) techniques. 

Figure^ shows the RDFy(r) functions (thin lines) ob- 
tained by Fourier transformation of the Sij{K) factors. 
From this picture one can see that the RDFNi-Ni {r) func- 
tion shows a very weak first coordination shell, located 
at about r = 2.64 A, while the second and third ones, 
located at about r — 4.5 A and 6.7 A, respectively, are 
very intense. The RDFi\fi-Zr (j^) function displays a well 
isolated first coordination shell, located at about r — 2.77 
A; a small shoulder at about r = 3.83 A, and well defined 
and intense second and third coordination shells, located 
at about r = 5.4 A and 7.9 A. The RDFzr-Zr(?') function 
shows the first coordination shell well isolated, located at 
about r — 3.3 A; the second one splits into two subshells, 
located at about r = 5.18 A and 6.35 A, which are very 
weak. The interatomic distances and coordination num- 
bers for the first neighbors obtained from these functions 
are listed in Table 
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FIG. 6: RDFNi-Ni(r), RDFNi-Zr(r) and RDFzr-Zr(r) functions 
obtained from the matrix inversion (thin lines) and RMC sim- 
ulations (thick lines) techniques. 



2. Partial Sij{K) obtained from the RMC simulations 

In order to perform the simulations we have consid- 
ered a cubic cell with 1800 atoms (600 Ni and 1200 
Zr), S = 0.01, and a mean atomic number density 
po — 0.054719 atoms/A'^. The minimum distances be- 
tween atoms were fixed at the beginning of the simula- 
tions at TNi-Ni = 2.2 A, TNi-Zr = 2.4 A and rzr-Zr — 2.6 
A. To make the simulatimis we used the RMC programs 
available at the InternetEj. The S{K, E) factors obtained 
at the photon energies 8330 eV, 17398 eV and 19200 eV 
were used as input data. These S{K, E) factors were 
cut at the same K value. The S{K, E) (solid lines) and 
S^'^'^{K,E) factors (square hnes) are shown in Fig. |^ 
and they show a good agreement. 

The S^j^^{K) factors (thick lines) are also shown in 
Fig. |5| together with those obtained from the matrix 
inversion method. The agreement among them is very 
good, mainly for the Szt-Zt{K) factor. For instance, the 
shoulder previously observed in iSzi-Zi (^) factor and the 
peaks located around K = 2.58 A~^, 4.3 A~^ and 6.2 
A~^ are also seen in the S^^^j^{K) factor; the minimum 
at if = 2.1 A~i and peaks around K = 2.7 A~\ 5.2 
A~^ and 7.6 A~^ observed in the 5Ni-Zr(i^) are also seen 
in the the 5^^zj(if) factor; and it is interesting to note 
that the minimum previously observed in the iSNi-Ni(if) 
factor is now replaced by a well defined maximum in the 
S^^^{K) factor. Other peaks observed before in the 
Sm-m{K) factor are also seen in the 5^j^(if) factor. 
The RDF^^*-^ (r) (thick lines) functions are also shown in 
Fig. 1^. The interatomic distances and coordination num- 
bers for the first neighbors extracted from the G^j^^ir) 
and RDFg"c:(-^) functions are also listed in Table ||. 

It is interesting to note that although we have used 
eight independent S{K,E) and four DSFi{K, Em, En) 
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1 2 3 4 5 6 7 8 



FIG. 7: S{K,E) (full lines) axiA S^^"^ [K , E) (squares) factors 
for the amorphous NiZr2 alloy. 



factors to obtain the three Sij{K) factors by the ma- 
trix inversion method, the iSNi-Ni (K) factor shows a min- 
imum between K — 1.5 A"^ and 3.5 which have no 
physical meaning. On the other hand, using only three 
S {K, E) factors as input data in the RMC method we 
have obtained excellent results, even for the ^Ni-NiC-f^^) 
factor. 1^ 1^ 

In a previous study, Lee et aln and Mizoguchi et ala 
reported the Sij{K) factors obtained for the amorphous 
NiasZrgs and Ni36Zr64 alloys, respectively. Lee et abused 
alloys containing mixtures of Co and Fe in small quanti- 
ties while Mizoguchi et al used the isotopic substitution 
method. A comparison among the Sf^^^{K)^ G^j^^{r) 
and RDF?!-'^^(r) functions obtained in this work, by us- 
ing only one sample, with those reported by them shows 
that they are in a good agreement. 

By defining the partial bond-angle distribution func- 
tions Qi-j-i{cos6) where j is the atom in the corner, 
we calculated the angular distribution of the bonds be- 
tween first neighbor atoms. The six Oi_j_;(cos0) func- 
tions found are shown in Fig. |^. This kind of information 
cannot be obtained by the matrix inversion method be- 
cause the partial distribution functions give only a one- 
dimensional description ofJJae atomic arrangement. 

Hausleitner and Hafnercil investigated several amor- 
phous alloys formed by transition metals using molec- 
ular dynamics simulations. They obtained the structure 
factors, coordination numbers, interatomic distances and 
also the bond-angle distribution functions. According to 
them, if the components of an alloy have a large difference 
in the number of d electrons there is a pronounced non- 
additivity of the pair interactions and a strong interaction 
between unlike atoms. Therefore, in certain amorphous 




cos 9 cos 9 



FIG. 8: Oi_j_;(cos^) functions obtained by the RMC simu- 
lations. 



alloys the interatomic distance of heteropolar .bonds can 
be shorter than the distance of homopolar onescJo. This 
also favors the formation of trigonal-prismatic units. The 
shortening effect was not observed in this study. The 
bond angle values of triangular and square faces of a 
distorted prism are found around 9 = 60°, 90°-100°, 
109° and 147°. Thus, the bond angle values shown 
in Fig. ^ suggest that distorted polyhedral units with 
triangular faces are present in the atomic structure of 
the amorphous NiZr2 alloy and if there are square-faced 
polyhedral units, they are found in a very small quan- 
tity. The presence of this kind_.of units was already vsi^ 
ported in anaqrphous Ni6oTi4oL3, NiasTigTc!!, Pd82Sii8c£l 
and NiesBgsB aUpys. 

Havinga et a/M3 reported crystallographic data for sev- 
eral compounds with CuAl2 (C16)-type structure. We 
have used their data for the NiZr2 compound in the 
CRYSTAL OFFICE98 software^ to obtain the inter- 
atomic distances and coordination numbers for the first 
neighbors listed in Table ||. Its structure has angles 
around 6 = 53-56°, 64°, 96-112° and 143-154°. For 
the amorphous NiZr2 alloy the 8i_zr-/(cos 6*) functions 
{i and I can be Ni or Zr atoms) found are very similar 
to those present in the NiZr2 compound. On the other 
hand, the Oi_Ni-;(cos0) functions show some important 
differences. The Ni-Ni-Ni sequence in the compound is 
linear [6 — 180°) while in amorphous NiZr2 alloy it can 
also be found at triangular angles {9 = 60° and 112°); 
the Zr-Ni-Ni and Zr-Ni-Zr sequences found in the amor- 
phous NiZr2 alloy show bond angles around 9 = 180° that 
are not found in the compound. These data suggest that 
in the amorphous alloy the local structure around a Ni 
atom can be more distorted than around a Zr atom. 

The chemical short-range order parameter is a quanti- 
tative measurement of the degree of chemical short-range 
order. There are several definitions of chemical short- 
range order parameter^. We have adopted the War- 
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ren chemical short-range order (CSRO) parameter ait^l 
given below to compare the amorphous NiZr2 alloy and 
the NiZr2 compound. 



ai = 1 — 



C2{ciN^+C2Nl) ' 



(5) 



Here, N^^ = A^Ni-Zr and = iVNi-Ni + ^Ni-Zr and 
N2 — Nzi-Zi + -^Zr-Ni are the total coordination numbers 
in the first shell. By using the coordination numbers 
given in Tablcp, we found af^^ = -0.029 and for the 
compound ctj^ ' — —0.024. The chemical short-range 
order in the amorphous NiZr2 alloy is similar to that 
found in a solid solution, for which ai = 0. The same 
behavior is observed in the NiZr2 compound. 



V. CONCLUSION 

The amorphous NiZr2 alloy was investigated by using 
AWAXS, DAS and RMC simulations techniques. The 
three Sij (K) factors were determined by using two in- 
dependent ways: a combination of AWAXS and DAS 
techniques and by combining AWAXS and RMC sim- 
ulations techniques. The iSNi-Zr(-^) and iSzr-Zr(-f^) fac- 
tors obtained by both combinations showed an excellent 
agreement, but the simulations gave the best result for 
the iSNi-Nil^i^) factor. The coordination numbers and in- 
teratomic distances calculated from the RDFij(r) func- 



tions obtained by these two combinations show a good 
agreement among themselves and with those extracted 
from the deconvolution of the DDF^i (r) and DDFzi- (r) 
functions. 

The results achieved in this study suggest that 
the combination of the RMC simulations method and 
AWAXS can be used to obtain more reliable Sij (K) fac- 
tors than those obtained by the matrix inversion method, 
in particular for the Ni-Ni pairs, which have the smallest 
weighting factor. This was already observed iii,a neutron 
diffraction study performed in a CuBr alloyc3, which is 
also an ill-conditioned problem similar to ours. In addi- 
tion, this combination gives us some evidence of the pres- 
ence of distorted polyhedral units with triangular faces 
in the atomic structure of the amorphous NiZr2 alloy. 

The calculated Warren CSRO parameter indicated 
that the chemical short-range order found in the amor- 
phous NiZr2 alloy is similar to that found in a solid so- 
lution, and there is some resemblance with that found in 
the NiZr2 compound. 
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